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CROUP-KINETIC THEORY OF TWO-PHASE TURBULENCE 


C.M. Tchen 

Department of Mechanical Engineering 
City College of New York 
New York, N.Y. 10031, USA 


Abstract 

The two phases are governed by two coupled systems 
of Navier-Stokes equations. The couplings are nonlinear. 
These equations describe the microdynamical state of tur- 
bulence, and are transformed into a master equation. By 
scaling, a kinetic hierarchy is generated in the form of 
groups, representing the spectral evolution, the diffusi- 
vity, and the relaxation. The loss of memory in formula- 
ting the relaxation yields the closure. The network of 
sub-distributions that participates in the relaxation is 
simulated by a self-consistent porous medium, so that the 
average effect on the diffusivity is to make it approach 
equilibrium. The kinetic equation of turbulence is derived. 
The method of moments reverts it to the continuum. The 
equation of spectral evolution is obtained, and the trans- 
port properties (eddy diffusivity and eddy viscosity) are 
calculated. In inertia turbulence, the Kolmogoroff law 
k 5/3 for the weak coupling and the spectrum k * 1 for the 
strong coupling are found. The density fluctuations have 
similar spectra. The numerical coefficients are determined 
analytically* 



1. Introduction 


The motion of small particles in an incompressible 
turbulent fluid is important for many applications (Soo, 
1967) , and from the fundamental point of view has at- 
tracted the attention of many authors as a statistical 
problem of many bodies. The earliest treatment of sus- 
pension (Tchen, 1947) used the Langevin equation (Langevin, 
1908) 

d - (t) + * v v(t) = X v (t) (1) 

dt — — 

for the motion of a single particle, where v(t) is the 
velocity fluctuation, a v is the coupling coefficient, and 
X y (t) is the driving force. A second Langevin equation 
may also be written for the Lagrangian representation of 
the modified fluid. 

The driving force is due to the drag that is experi- 
enced by the particle from the moving fluid. It is pro- 
portional to the fluid velocity u(t) as a v u(t), and may 
include a Basset memory. The interaction of many particles 
calls for a system of many coupled Langevin equations 
(Herczyfiski and Piefikowska, 1980; Mazur, 1982; Deutch 
and Oppenheim, 1971; Ermak and McCammon, 1978; Ramshaw, 
1979) . 

From the Langevin equation with a drag and a memory, 
Tchen (1947) calculated the Lagrangian correlations of 
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velocities and found the equality 



v(t)v(t-T: )> 



A fluctuating function 


(2) 


^ — V 

()=() + 0 


(3) 


can be decomposed into an ensemble average 

0 s < > (4) 

and a fluctuation (~). The operators 

A « A + A (5) 

may be used for the scaling, where A=1 is the unit opera- 
tor. 

By the use of the evolution operators 

U v (t,t-T ) and U u (t,t-r) (6) 

to specify the trajectories of the v-fluctuation and the 
u- fluctuation, respectively, we can write the Langevin 
equation (1) and the equality (2) as 
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(7) 


^ M /VA 'V 

- — + oc v(t) = AU ( t f t- X ) X ( t— T ) 

dt ~ - 

and 

) 

dT^(t,x)AU v (t,t-X)v(t-X))> = /dr^u(t,x)AU (t^t-T J'utt-T )) (8) 

The fluctuating operation A makes the quantity which 
follows to be a fluctuation. The force X v (t-t) in (7) 
is registered at the time t - t and at the position X(t-t) 
by an observer riding on the particle that follows the 
perturbed trajectory of the . v-motion, while at time t, 
the force is 


X (t) = X (t ,x) . (9) 

v « 


The integrals 

J d-o (v(t,x)AU (t,t-T )v(t-r )) 

0 ~ N. V — 


( 10 ) 


and 


I, 


Of 

dr 


<^u(t,x)AU u (t,t-T )u(t-t ))> 


(ID 


can be called dif fusivities. The correlation functions are 
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I 


calculated from the velocity fluctuations at the two in- 
stants along their own trajectories, as specified by the 
evolution operators (6) . 

With this distinction between 


U v (t,t-T )u(t--C ) and U u (t,t-T)u(t-T ) , (12) 


the quantity 


i. 


CD 

< M -V A 

U (t ,x) AU^ 


(t,t- r) 


u(t-r)') 


(13) 


as formed by the u-f luctuations along the v-tra jectory , 

i "■ — 

is not a diffusivity as a property, but is a functional 
of the diffusivity (11). In this sense, the Langevin 
| equation (7) and the kinematic relation (8) are in fact 

nonlinear. The nonlinearity arises from the transforma- 

i 

tion between the evolution operators (6) . 


The diffusivity in the operator form calls for the 

I 

need for a Lagrangian-Eulerian transformation (Pismen and 
Nir, 1978; Reeks and McKee, 1984; Gitterman and Steinberg, 
1980; Gouesbet et al., 1984). It soon becomes clear that 
without the dynamical knowledge of turbulence, the trans- 
formations between the Lagrangian and Eulerian representa- 
tions and between the two evolution operators can not be 
successful . 


1-5 


The treatment of many interacting particles and the 
transport theory of the stochastic system belong in a 
natural way to the many-body statistical mechanics (Bogo- 
liubov, 1962). Yanko (1980) has suggested Bogoliubov's 
method for treating the suspension of particles. The 
method is best suited for a system of large numbers of particles 
that interact through their own potential of interaction, 
e.g. the Coulomb potential in plasmas, but is not well- 
developed for the interaction among the particles in a 
flowing fluid. It would be more pertinent to describe 
the interaction by a system of many Langevin equations, 
for the many particles embedded in the incompressible 
flow. 

If we are not concerned with the microdynamics of 
the interactions and the microhydrodynamics of the in- 
trinsic viscosity in our theory of turbulence, we can 
choose an effective friction constant, such as the Darcy 
damping for the ensemble of Langevin equations treated 
as a porous medium. Through the transformation into the 
Fokker-Planck equation, and the subsequent reduction in- 
to the continuum, we will obtain two systems of Navier- 
Stokes equations for the two phases, with the effective 
coupling constants as given parameters. These systems 
describe the microdynamical state of two-phase turbulence. 

For the development of a theory of turbulence, the 
direct hydrodynamic method has shown the difficulties of 
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the closure of hierarchy from the start of the problem 
(Genchev and Karpuzov, 1980; Elgobashi and Abou-Arab, 

1981). The statistical method by Bogoliubov's approach 
with the closure at the triplet distribution seems un- 
satisfactory for turbulence. For this reason, we develop 
a group-kinetic method to derive the kinetic equation for 
the distribution function in the macro-group. The closure 
is found by memory loss. We investigate the eddy viscosity 
and analyze the spectral structure of turbulence. 

It is generally believed that a small amount of 
particles present in the fluid will drastically reduce 
the drag of the fluid. The laboratory measurement of the 
turbulent intensity indeed shows a decrease of fluid tur- 
bulence with the increasing concentration of particles 
(Elgobashi and Abou-Arab, 1983) . Other measurements of 
the energy spectrum in the inertia turbulence confirm the 
Kolmogorof f law k‘ 5/3 , but show unexpectedly an increased 
spectral level (Goldschmidt, 1972) . To lift this paradox, 
we investigate the inertia turbulence with strong coupling 
and weak coupling. In the direction of increasing wave 
numbers, our results show that the spectral law k 1 for 
• strong coupling precedes the spectral tail of power 
law k~ 5/3 for weak coupling. The power law k* 1 is reduced 
by the particles and the power law k‘ 5/3 is increased by 
the particles. The tendency of this shift from reduction 
to increase is indicated by the lidar measurement in the 
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atmosphere before and after the rain (Gurvich and Pokasov, 
1972) . 


2. Fundamental equations of two-phase flow 

We consider two fluids, called the two-phases v and 
u. Both fluids are governed by the Navier-Stokes equations 
with coupling, in the form 

f v ( \+v*V ) v - - + y v n v (u - v) 

\ i + '7- f v ? “ 0 ' 

for the velocity v, the pressure p v , and the density p v , 
and 


(14) 

(15) 



u* V ) u 



(16) 


* 

t 



+ V 



= o 


(17) 


for the velocity u, the pressure p u , and the density p u . 
The coupling constants are Y v /Y u * 


By definition, the mass densities 


?v ■ m v n v ' f« “ m u n u 


(18) 


can be written as the product of masses m v , m u , with the 
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number densities n v » n u . The kinematic viscosities v v , v u 
are not written here since they are negligible as com- 
pared with the eddy transport properties, but will be in- 
cluded in the energy considerations. 

The two fluids are incompressible, but the coupling 
depends on the density. In investigations of this kind, 
we allow for the variation of density insofar as it modi- 
fies the coupling only. This procedure is called the 
Boussinesq approximation (Chandrasekhar, 1983) . Thus, in 
! eqs. (14) - (17) we replace p„, p„ by their stationary values 

i V U 

P v / P u and reduce the hydrodynamic systems into the follow 
ing: 

(i) for the v-phase 

( ^ t+ v. F +J v )0 = X v (19) 

V- v = 0 (20) 

— -w 

I 

< V l' 7 >"v 3 V N,v 7 \ ' < 21 » 

with 



A 



* A 

+ U 
V 



( 22 ) 
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(ii) for the u-phase 

(> t + u-7+« )u = x 

7-G = o 

<^t - U’3>“u - V N,u^u 

with 


(23) 

(24) 

(25) 


A ^ A ^ 

X = E + v v , E 
u u u ' u 


1 ^ 

— 7 P 

p *v 
/u ^ 


(26) 


The coupling coefficients 


<* = N , o< = v N 
v v v ' u u u 


(27) 


fluctuate by the fluctuating number-densities, and the 
stationary coupling coefficients are 


c< = "Y /m , p< = V /m 
v 'v 7 v u u' u 


(28) 


The normalized densities 


A .* 


N =n / n , 
v \r v ' 


A vA 

N =n /n , 
u u 7 u 


(29) 
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are the densities n v ,n u normalized by the stationary values 
n v , n u . The stationary quantities are written without the 
symbols (3). 

It is to be remarked that the Navier-Stokes equations, 
in the form (14) -(17) or (19) -(26), are valid for two gases 
of number-densities (29) or for two liquids of densities 
(18) . In the special case of the suspension of particles 
in a liquid, where the coupling is due to the friction of 
the particles only, we have the simplification 

N u - N v , (30) 


without altering the form of the basic equations (19) -(26). 

In view of the symmetry, it suffices to treat the 
first system and omit the subscript v, unless needed for 
clarity. The results will equally apply to the second 
system by interchanging v and u. 


3. Group-kinetic method 
3.1. Master equation 

The inhomogeneous partial differential equation (14) 


can be transformed into the homogeneous partial differen- 
tial equation 
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called the master equation for the distribution function 
f(t,x,v), with the differential operator 


A 

L 


V. v 


+ 


A 


+ 



(32) 


and the notations 


= "^ / *> t / ^ = V ^ v 


(33) 


For establishing the self-consistency between the 
master equation and the hydrodynamics of the two phases 
under the Boussinesq approximation, we write 


f ( t , x , v) 



- v(t, 


?] 


(34) 


and 


f (t,x,v) 


N(t,x)^ fv-v(t,x)] 


(35) 


for (19) and (21), respectively. 


3.2 Group-scaling 

The master equation describes the microdynamical state 
of turbulence in too many minute details. A smoothing pro- 
cedure is necessary for a statistical treatment. By the 
Reynolds decomposition (3) , we divide 

A — ^ 

f - f + f (36) 

into an ensemble average f and a fluctuation f. The group- 
scaling by 


0 

A = A + A’ 

A' = A 1 + A" 

decomposes the fluctuation 

0 

f = f + f' 

into a macro-group f° and a micro-group 
f' = f 1 + f" 


(37) 


(38) 


(39) 


(40) 


The micro-group is re-decomposed into the first-order 
group f 1 and the group 
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f 


f (2) + f (3) 


• • • 


(41) 


+ 


of sub-distributions. 

The three groups 

f # , f' , f" (42) 


represent the three transport processes of evolution, trans- 
port property, and relaxation. The evolution of f° occurs 
in a medium that offers a transport property as shaped by 
the micro-fluctuation f' . Subsequently, the transport 
property approaches its equilibrium (t -*■“») by a relaxation 
from f" . The relaxation is a complicated process obtained 
by following the perturbed trajectory and by closing the 
memory chain. The distinction of the transport processes 
provides a physical demarcation in the groups. Mathe- 
matically, the groups are sets of Fourier components, so 
that the intensities of the field fluctuations 


rk f 00 f* 

(x' 2 )=2j dk'Gfk') ,(x'^=2£ dk"G(k") ,<X ,,2 >* ^dk^GOc'") 


( 43 ) 


are portions of the spectral distribution G(k) within the 
wave number intervals 


(0,k), (k,<*> ) , and (k",®), 


( 44 ) 


respectively. 
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The three groups of increasing orders by the scalings 


A , A' , A" 


(45) 


do not refer to decreasing magnitudes, but to decreasing 
coherence. Thus, the decreasing times of correlation, by 
the inequalities 




(46) 


will introduce the quasi-homogeneity and quasi-stationarity 
at successive group levels to limit the interactions to 
nearest-neighbour s . 

3.3 Hierarchy and closure 

By A 0 and A 1 , we scale the master equation (31) into 
the following equations for the evolution of groups: 


O o 

(\+AL)f = 


0 — *11 
- L f - A Ift 


(47) 


(V^f 1 -LW) >AV 2 ) f (2) , ... 


(48) 
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and the sequence continues in a hierarchy. The 
sub-distributions form a global network (41) like a 
lattice. 

On the left hand sides of (47) and (48) , the 
advections can be written as 


o * 0 0 o*o 

ALf = Lf-ALf (49) 

0 


1 T l f l 

A L f = 


f 1 - 


A 1 L 1 


(50) 


The notations 

L o s L + L , L x = L + L° + L 1 (51) 

are used. By the approximation of nearest-neighbour 
interactions, the last terms 

- A* L 1 f 1 = - A*L' f* = £' j f ° } (52) 
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-AV 2 >f< 2 > 


(53) 


in (47) and (48) are called collision and sub-collision, 
respectively. Since 


S'f } ' 


are operators in v, the memory continues. 

Physically, the evolution of f° proceeds in a col- 
lisional environment where the fluctuations by f' can 
shape an eddy viscosity. For its approach to equilibrium, 
this transport property will need a relaxation, where we 
must find a way to cut the memory and thereby to close 
the hierarchy. 

By an inspection of the chain of collisions in (52) 
and (53) and of the linkage among the sub-distributions 
in the network (41) , we visualize that the memory is at 
the same time essential for maintaining a non-Markovian 
behaviour in the kinetic equation of turbulence, and 
should ultimately vanish for obtaining the closure and 
the irreversibility. It can be seen that the memory con- 
tinues through the v-dependence among the sub-distribu- 
tions linked in the network (41) . In fact, it is the 
averaged statistical effects that play the role of relaxa- 
tion and controls the evolution of f l . Consequently, we 
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can replace the network by a porous medium that offers a 
friction Q " , called the Darcy damping (Brinkman, 1947). 
This amounts to replacing £"{ } in the role of operator 
by g" as a function independent of v, i.e. 

S'{) a 5" . (55. 


This function will be determined self-consistently . In so 
doing, the sub-distributions in the network lose their 
individuality as v-dependent and thereby are not capable 
of linking the memory in the homogeneized lattice that 
turns into a porous medium. 

4. Kinetic equation of turbulence 

By substituting for the sub-collision (55) , we trans- 
form (48) into the form 


( ^A 1 !) f 1 = ” I^tf+f*) +S" f 1 


(56) 


or 


= - L 1 (f+f d ) , 


with the differential operators 


A ✓N v p A, II 

= v V + <K + X- 6 and JL- L - 


(57) 


(58) 
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The two operators differ by the Darcy damping for 
representing the porous medium at large. 

It is to be remarked that the inhomogeneous equation 
(57) can be integrated by means of a kernal t^(t, t - t) , 
called the evolution operator, that satisfies the homo- 
geneous equation 


fjL + l (tj) $(t,t x ) = 0 , 


(59) 


"D t 


with the condition T*£(t,t) =1. The solution is 


li(t 


,t x ) = exp [-/ dt' x £ (t L )] . 

t. 


(60) 


By means of this operator, we integrate (57) to get 
the distribution 


A f - 0 

f 1 = -] dr A 1 U(t,t-T jJl 1 (f+f° )] , 


(61) 


and calculate the collision (52) as 


1.1 


tfff } = - A Lr 


= A^dT L 1 (t,x , v) A L (1( t , t-I ^L 1 (f+f*) | 


(62) 


A final substitution into (47) yields the kinetic equation: 
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( 63 ) 


0 0 — . / ^ \ 

( W f = " Lf + 5 j f <*■£>} • 

This equation differs from the Bogoliubov equation and 
its generalization to include the memory (Tchen, 1959) , 
and will serve as a basis for deriving the transport 
properties and the spectral structure without the need 
of the accompanying pair-distribution function. 

It is worth noting the following properties of the 
collision {f 0 }: 

(i) As a collision operator £' { } is deterministic 
and stationary, so that the upper limit of integration in 
(62) can be replaced by 00 without altering the value of 
the integral, by 


T v > T v , (64) 


from (46) . 

(ii) The nearest-neighbour interaction approximation 
between groups is valid by (46) , so that the correlation 
of L 1 in (62) can be replaced by a correlation of L' . 

(iii) The approach to equilibrium depends on the loss 
of memory by relaxation, so that we can write 

A __ 

U (t,t--T) = U(t,t-T) . (65) 
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By these approximations, we transform the collision 


(62) into 



where the scaling gives 

• - « « 

A (f + f ) - f . 

By the commutability 

X' (t,x) » > *x' ( t,x) , 

- *- - — 

we can write the correlation 

(v (t,x,v)*l?(t,t-T)L' (t-Z)> = ^ • /x' (t,x)5(t,fc-T)X' (t- 

+ x 2 <^' (t,x)ii (t,t-T)N* (t-t)) 
and the collision operators 

c'O-ciU* «n 

in two components 

— 


/ ( 66 ) 


(67) 


( 68 ) 


(69) 


(70) 


(71a) 
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and 



9 


with the diffusivities 






(t,x)U(t,t-T)X' (t-^ 



(71b) 


(72) 


(73) 


It is to be remarked that the auto-correlation of 


L' 


X N' + X' 


.'J 


(74) 


will entail the auto-correlation of X' and of N', while 
the cross-correlations are negligible. 

In conclusion, the kinetic equation (63) with the 
collision ( 55 ) will serve as a basis of a transport 
theory. 


5. Hydrodynamic equations of turbulence and cascade transfer 
5.1 Kinetic form of cascade 

The hydrodynamic equations (19) and (20) describe the 
microdynamic state of turbulence. They were transformed 
into a master equation as the starting point of our kinetic 
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theory. By group-scaling and the closure by memory loss, 
we derived the kinetic equation of turbulence (63) . By 
the moment method, we can convert the kinetic representa- 
tion back into the continuum to find the hydrodynamic 
equations of turbulence, as follows: 


(• j t + v ‘- 7-* 2 


0 

+ A<)v 



o 

+ J 


(75) 


and 


o 

"• V - 0 . ( 76 ) 

It is to be stipulated that the macro-kinetic equation for 
f° is an irreversible equation, and does not allow to be 

A A 

reverted from f into f to generate higher-order groups 

A 

of moment, as did the reversible master equation for f. 

The driving force is 

o 9 0s> 9 o i 0 

X = E + A <* u , with E = - — V p , (77) 

- - _ y - 

and the hydrodynamic force due to the collision is 

J° = [ dv v SM f°\ • (78) 

J — “ [ J 

We have added the kinematic viscosity v. 
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By (70) -(73), we can write 


J = J x + J, 

w» 


( 79 ) 


with 



(80) 


(81) 


Upon multiplying by v° and averaging, we get the 

•N 

energy equation 

i \<:f 2 > = - t " + «' - 1* -z" ■ (82) 

We have assumed homogeneity and isotropy, so that 




V •< 


o o\ / « 
v p / - \ p 7 




= 0 


* 


(83) 
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The production function 


W 



( 84 ) 


and the damping function 

are originated from the couplings. The dissipation 
function is 

^ - V ((7vV) . 


( 85 ) 


( 86 ) 


Finally the transfer function 


t 


T 


= - <V • j 4 /* 

= - v'.fd vv^'j. 


( 87 ) 


It is the result of interaction between the macro-dustribution 
and the micro-field fluctuations. It governs the cascade 
flow across the spectrum. 

The inertia turbulence is governed by 


= l , 


(fit) 
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where the energy dissipation 



(89) 

is increased by the drag a. 


From the collision of two components (71) , 

it entails 

the transfer function 


0 * o 

T - T + T 

N A X 

. (90) 

of two components 


t n - - ■ z 6 n O 


- - » <v‘ 2 > 

(91) 

and 


T x = - ] d z( v ° • ; 5 xf f '} 


= - f dv ^ v • vt)*^ x •bjf 0 j y 

(92) 


5.2 Fluid form of cascade 

From the scaling of the fluid equation (19) by A 0 
and h' , we get the equation of evolution for the macro- 
velocity 
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*__v » ^ O v O A 0 

( ^ +v *V)v =X - A X v + J 


( 93 ) 


The force of dissipation 


1 m J X + J N 


( 94 ) 


has the two components 


o o 

J v = - V • A v* v' 


( 95 ) 


J N = -ei A N' V' , 


( 96 ) 


in contrast to the kinetic form given by (80) and (81) . 
These are due to the stresses 


-A v' v' 


(97) 


and 


-A N' v* 


( 98 ) 


from the micro-fluctuations. They represent the statistical 
effects of the micro-fluctuations upon the evolution of 
the macro-momentum. 

Upon multiplying (93) by v° and averaging, we get the 
energy equation in the same form as (82) , with the transfer 
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function also in two components now in the fluid form: 


f ® V 

T x =<v *7* A v* v'/ (99) 

T N *^\ v * A n ' v '/ / (100) 

in contrast to the kinetic form given by (91) and (92) . 


5.3 Transport of the gradient type and transport of the non- 
gradient type 

The flux (97) 


o 

- A v' v' 

j i 


K! V u. 

JS SI 


( 101 ) 


for the transport of v[ by the velocity fluctuation v^ is 
proportional to the velocity gradient V s u?. The transport 
coefficient Kj s is called eddy viscosity. This transport 
relation was used by Heisenberg (1948) to represent the 
force of dissipation 


J* = K ' vV (102) 

in the notation (95) W ith9 rou P s * It served to determine the 
transfer function 


o 



K' ^(7uV^). 


(103) 
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The same relation for this transport of the gradient type 
is called the Onsager relation in the thermodynamics of 
irreversible processes. These phenomenological arguments 
can not lead to the determination of the transport co- 
efficient, for which a kinetic approach is necessary. 

This is true for the eddy viscosity as well as for the 
molecular viscosity. 

By comparing (102) with (80), or (103) with (92), we 
find the relationship in the form 

K' 7 2 (104a) 


and 


^=K"^ 2 , (104b) 

between the collision coefficient and the transport co- 
efficient, when they are independent of v, i.e. after the 
memory cutoff. Therefore, this relationship identifies the 
collisions in the calculation of the relaxation, and en- 
ables the self-consistent determination of the eddy vis- 
cosity. 

On the other hand, the transport 

A N' v' = o< v , (105) 
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unlike -A°v'N’ , is not of the gradient type, as seen by 
comparing (96) with (81) . 


6. Perturbations of the trajectory 
The Lagrangian correlation 


<^X' (t,x)li (t,t-T)X' (t-T)^> = /x' (t,x)X’[t-r x(t-T)]^ (106) 

is defined as the correlation between two fluctuations at 
two instants (t,t-x) along the trajectory. The variable 
position 


x(t-r) 


= x - vX 


- J (z! 


consists of an initial position 


x(t) = x , 


a free flight vt , and a perturbed path 


/ (r) = i (r) + i (r) . 


(107) 


(108) 


(109) 
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By Fourier decomposition, we have the fields 


X'(t,x) - ffdto ' dk' X'(<i',k') e -i(^’t-k'*x ) 
X' [t-T, x (t-T)] = JJ dti” dk" X’ («",k")e~ 1<J " 


(t-f) ik"-x(t-t) 


and the correlation 


( 110 ) 

( 111 ) 


<^X' (t,x)X‘ [t-r,x(t-I)]^ <4* dk'JJW dk"^X , (g ' , k ' ) X 1 (<o " ,k" )^ 

e -i[( U 'W)t - (k'+k")- S J h0 . (k „ (VK (U2) 


The orbit function is 


h(T,k" /V ) 


A 

_ -ik"-vT /■ -ik"«J? (r)\ 
- 6 ^ ^ \^6 ^ « / 


(113) 


by (111) and (107). 

If the Fourier decomposition is made in an interval 
of time 2 T and an interval of length 2 L in three dimen- 
sions within which the field is stationary and homogeneous, 
we can average (112) in time and space, integrate it with 
respect to t, to obtain the diffusivity 

f drffdg" dk" (<■>", k")X’ (-o",-k"£> h(T,k",v), (114) 

^ X J 0 ^ " -* ~ ~ 

as defined by (72) . The factor of Fourier truncation is 
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y = 1 / 3 ) 3 

X TL T \ Ly . 


( 115 ) 


With a free flight in a frictionless medium, the 
dif fusivity 




d«J" dk" k")X’ (-«; k")) l (<j"-k".v) 


has a resonance, as is common in linear instability. In 
weak turbulence, the trajectory of X' is perturbed by a 
smooth streaming from 1° -fluctuations. 

We deal with a strong turbulence by neglecting ui" 
in (114) and reduce it into the form: 

= f tJ> d? fdk" (k")X' (-k")\h(T,k",v) , 

X J o J ~ ' - - “ - 7 - - 

where x = ( tt/L) 3 is the factor of Fourier truncation with- 
in the interval of length 2 L in three dimensions. 

The orbit function (113) has a component 
( e -ik”.|(r)> . J d j .W p(T J, 

from the path fluctuation Mt), and can be calculated by 
means of the probability of transition p(x,i). A theory 
of transition can be developed on the basis of a master 


(116) 


(117) 


(118) 
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equation by a method similar to the one presented earlier. 
As a result, we will equally have a Darcy damping to 

represent the friction by the porous medium: 


-ik"* >(r) 


!> 



(119) 


Analogously, we have 

-ik"* Jt* (x) 

with ^ 1 = , by the approximation of the nearest-neigh- 

bour interaction. Finally, the coupling constant a enters 
as an external parameter, giving 


= e 






A 


( 120 ) 


'k 

u 


e 


( 121 ) 


By collecting (119) -(121), we find the orbit function 


h * h * h v h D h Jt ' 


( 122 ) 


with the following components: 


h v (-r,k",v) = 


-ik" • v T 
0 *•* -* 


(123) 
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h D (x,k") 


(124) 


h^(-c,k") - e . (125) 

By (104) , the last two orbit components can be written in 
the form: 


h D (T,k") 




exp (-k K"T ) 


(126) 


h (T,k”) 


exp (-k K'T 

Z 


) 


(127) 


Note that 1 = 0, and that i° (t) is neglected, because 

s 

the evolution operator 


A' U (t,t-T) (128) 

in (61) selects the micro-group by A'. It may play a role 
of streaming 

The only component that represents the path-perturba- 
tion is h 4 (T,k"). It performs the physical role of relaxa- 
tion for the eddy viscosity K' as a transport property to 
approach the equilibrium at a self-consistent rate k" 2 K‘. 

I 
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7. Cascade transfer by eddy viscosity 

By (92) and (71) , the transfer function t£ can be 
written in the kinetic form, as follows: 


T = 
A X 


- f dv v* (t,x) 5 jf }>} 

= ~J dv v«^v°(t / x) f (t“ T )J^ 

* - Jdv v *^ x j^ y (T?»v)j> 

= -^ x {Jdv v \ 2 y (T,V) | , 


(129) 


with 


» / « * \ 
'P (T,V) = \V (tx)f (t“T) / 


(130) 


and 


- tr S' 
* X 


In Fourier form, we write 


(131) 


(T,v) “J'dk 1 / (v° (k *) f *(-k' ,v)^> h(r,k',v) 


and 


(132) 
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*( w 


(133) 


I 2 *f (r,v) = fdk' <^v°(k') t (-k • ,v)^} 2 h (v,k • ,v) 

**» " J - - - •** 


The differentiation applies to the fast varying function 
h(T,k' ,v) and not to the slowly varying function f°(-k' ,v) 
as an approximation. In the same manner, we write 



(134) 


from (117) and (128). 

By collecting (133) and (134) , we transform (129) 

into 


| dcj<5k"%(x' (k")X* (-k")^> h(v,k",v*=0) 

* J'^Z. f ("k'v^Atr/k* ,v) 


Here we have two orbit components 


(135) 


h(T,k' , ,v=oj^ 2 h(-r,k’,v) = -k^T^frA'SvN)) hfr,k’,v) , 


(136) 


The second component can also drop its v- 

**• 

dependence after its differentiation by 3 2 , so that 

h( T ,k" ,v=0) h(T,k',v) = h(r,k\v=0) h(r,k',v=0). (137) 
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In the competition for the dampings, the second 

component is ineffective, and can be approximated by 


h(r,k',v=0) = 1 . 


(138) 


By the use of ( 136) - (138) , we transform (135) into 


dr xr h(T,k",v=0) 


=£ f dk"2( ^X» (k")X' (-k "))f 
J * ** ^ 0 

X J k' 2 /vV).vVk')) 


= K' R . 


(139) 


By separating the integrals into two parts, we recog- 
nize the vorticity function as 


R =Jdk' k ' 2 ^ v ("k'j^ , 


and identify the eddy viscosity as 


(140) 


K 


' =U jdk" / (x' (k") X' He" j) J < 


dt x h(r,k",v=0) 
0 - “ 


(141) 


or as 


K ' = K\' 2 \ ■ 


(142) 
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Note that the diffusivity as a transport property 


S' = (k")X‘ (-k"))J dr h(T,k",v=0), (143) 

w * 

-v 

by (117) , can be used as operating on t 2 . 


8. Eddy viscosity 

We have derived the formulae for the transfer func- 
tion T® , the eddy viscosity K', and the orbit function 
h(T,k",v = 0), in (139), (141), and (122), respectively. 

We show that the spectral evolution is characterized by 
the eddy viscosity K'. This transport property approaches 
its equilibrium by a relaxation process. The analysis of 
the perturbations of the trajectory results in the orbit 

function h ( x , k" , v = 0 ) . 

<*>* ^ 

From (141) we can distinguish two eddy viscosities: 


K' = tr 






j dk" x/x’ (k")X' He" 




dk" 



(k") 


(144) 


(145) 


1-38 


They control the cascade transfer and determine the relaxa 
tion, respectively. They are energized by the field fluc- 
tuations of spectrum G(k M ) , and are characterized by the 
two time moments: 



h(c,k",v=0) 


h(r,k") . 


(146) 


(147) 


The latter is the moment of h^ in its role of relaxation, 
in the form 



(148) 


so that, by a substitution into (145) , yields 


K' 



G(k") (k' ,2 K^ )" 3 . 


(149) 


The eddy viscosity K£, as appearing in both sides of (149), 
is equally stationary and homogeneous since it is a func- 
tion of k. By grouping them together, we simplify (149 ) 
into the form 



(k") 



k"' 6 G (k "' ) j 


1/4 


(150) 
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For the calculation of K', we distinguish between 
the weak coupling and the strong coupling in the follow- 
ing. 


8.1 Weak coupling a« k" 2 K^ 

In weak coupling we neglect a in (122 ) an< ^ (146) to 

get 

^( k") =J° dz x 2 expj-k' ,2 fK"(k")+Kj (k")]} 


= 2 jk” 2 [K’’(k")+K^(k")] } 

= 2 [ 2k'' 2 K' (k”) ] 3 • 


-3 


(151) 


We note that by its higher order in the group, K"(k") is 
qualified to play the same role of relaxation as does 
K^(k"). Use of (126) and (127) has been made. 

It is not difficult to demonstrate the relation 



by differentiating the two eddy viscosities (144) and 
(150) and by comparing them. 
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8.2 Strong coupling k" 2 K^ « a 

The calculation of G(k"), as defined by (146), is 
more involved. The function has two time scales 


oT 1 and [ k" 2 (K' +K") ] ” 1 


(153) 


to distribute. As an optimization for fast approach to 
equilibrium, we write 


Cj (k“) =/^ dr }& z ' expj-k" 2 |K" (k")+K£ (k”)]T'| , (154) 


or 



[ 


K"(k")+K‘ 

X 




(155) 


when K" joins with for the relaxation process by its 
qualification of belonging to a high-order group. The 
upper limit t of integration in (154) can be put to 00 
in view of the' condition 


* » 



(156) 


Upon substituting (155) into (144), we find 
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K* = 


(157) 


r® 


- 

3* y k 


dk" k 1 


„-4 




The proportionality a" 1 fulfills the Einstein condition 
in Brownian motion. 


9. Equation of state relating pressure and velocity fluctuations 

The transfer functions T® and T® as obtained in (91) 

N A 

and (92) depend on N' and X', respectively. For T®. , a 

X 

relation between X' and v' is needed. It is a kinematic 
relation without involving dynamical processes and is 
called the equation of state. Since T® describes a 

A 

momentum transfer at the exclusion of the density fluc- 
tuations, we can use the hydrodynamic equation (19) with 
constant coupling, i.e. 6t = a, as our basis of derivation. 
By taking the divergence, we find the Poisson equation 

= r , (158) 


with the source 


C* ** 

r =77: Aw 


= W : (v v - Cv Zy j 


(159) 


By Fourier decomposition, we transform the equation (158) 
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for the E-field into the following equation for the inten- 
sities, as 

(jm-EHo) - k’ 2 </r(k) r(-k)). (160) 

It is evident by (159) that the right-hand side will have 
quadruple moments. These will be decoupled into products 
of quadratic ones by the assumption of quasi-normality. 
This is legitimate, since the equation of state deals 
with a stationary, homogeneous, and isotropic state of 
fluctuations, and notwith a dynamic process. We find the 


pressure-velocity relation in the form: 


The details of the derivation which uses the Fourier 
method are omitted here. 

For the sake of brevity, it will be convenient to 
introduce the notation 


r* = fdk' k ,2 /(w # (k')- w°(-k'£> 


with w= (v,u,£,X) . We find 




< 


(162) 


(163) 
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and 


0 o 5 0 

«X - % + “v \ 


(164) 


from (160) and (22) , respectively. 


By substituting (161) into (163) and subsequently 
into (164) , we transform it into 


I <<> 2 + ,/ v 2 < 


(165) 


By differentiating with respect to k, we obtain the 
following relation between the spectra: 


G0c)=|r°F v 0c) + « v 2 F u (k) 


(166) 


This relation will serve in transforming the transport 
coefficients ' and K' that were written in terms of X 

X vt 

into those in terms of velocities. 


The functions G(k), F v (k) and F u (k) are spectral 
densities, giving the spectral intensities: 


Jdk" (k") *X* (-k"))= 2 J dk" G(k") 
J^dk'^^v°(k') *v°(-k'^= 2 J* dk' F v (k’) 

f dk'^^u 0 (k') *u* (-k'^= 2 J ^dk' F u (k') . 


(167) 


(168) 


(169) 


10. Spectral structure of velocity fluctuations 

10.1 Inertia turbulence with weak coupling 

The inertia turbulence is governed by the cascade 
flow at the constant rate of dissipation. If the coupling 
is negligible, i.e. 

V « k 2 ** , (R°) 1/2 , (170) 

the two phases are separated. 

The spectral flow is governed by 

T° = K’ R* = £ , (171) 

and is independent of a. The small-scale transfer is 
chosen, because we are treating the range near the dissi- 
pation end of the spectrum. 

The expression (150) for the eddy viscosity in 
terms of the spectrum G(k) of the field spectrum can be 
converted in terms of the velocity spectrum F v (k) by the 
equation of state (166) which is 

G(k) =|r* (k) F v (k) < 172) 

in weak coupling. We find 
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The numerical coefficients are predicted analytically, 
and are in good agreement with the best experimental 
values. 

The -5/3 law is identical to the Kolmogoroff law, 
except the rates of dissipation e y , e u are augmented by 
the drag coefficients, by (89) . Such an increase has been found in 
experiments, where the coupling coefficients are small 
to satisfy the condition (179) , but not so small as to 
be negligible in (89) (Goldschmidt, 1972) . 
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I 

10.2 Inertia turbulence with strong coupling 

In strong coupling, i.e. 

* >> , (rV /2 , (178) 

the spectral flow is governed by equation (171) , rewritten 
f as 



(179) 


for the v-phase. 


The vorticity function is 




dk' k' 2 F v (k') 


(180) 


The eddy viscosity is given by (157) in terms of the 
spectrum G y (k). It can be converted in terms of the spec- 
trum F u (k) by means of the equation of state (166), which 
becomes 


G v (k) = « 2 F u (k) 


in strong coupling. We get 


2 < v /~ 

K v = _ T J dk" k"~ 4 F u (k") [2K’(k")J _2 , (182) 

k ^ ^ 


1-47 



with 


K '/ = <3k " kM_6 *V (k " ) p v fc H >J V4 (183) 

by (173). 

The system (179), (180), (182) (183) governs the 

v-cascade. A parallel system governs the u-cascade, and 
is obtained by interchanging v and u. The two systems 
are coupled. 

By omitting the details of calculation, we obtain 
the following results: 


F (k) = A u 2 k 1 
v' * £ 

F (k) = A v 2 k' 1 

U * 

K - ' l W v < > 2 k ' 2 

with the friction velocities 


v t 3 ( V<v> 


1/2 


u 

i 


= < V*u> 


1/2 


(184) 

(185) 

(186) 

(187) 

(188) 


The numerical coefficient is predicted as: 
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( 189 ) 


A = 8/0 . 

«< 


The solutions can be obtained by a method of verification 
using a power law spectral distribution. 

1 1 . Spectral structure of density fluctuations 

The kinetic equation (63) is valid for density fluc- 
tuations, by definition (35). By taking the zeroth moment, 
we find the equation of evolution 


+ v°*V)n°= - v°-7N + 7 V+ j&XJ \ f ° (t_T) } 


(190) 


Upon multiplying it by N° and averaging, we derive the 
equation of spectral flow 


I \( N ° 2 /" *^ N ° V N - £ N - 0° . 


This is governed by the production function 


(191) 


-(n*v>7n. , 


(192) 


the dissipation function 


(193) 


C -*»<<*■'> 2 > 
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by molecular diffusivity v N , and the transfer function 



The transfer function 


0 =% 



has again two components 


©x ■ -J <T x {**Mr> 


( 194 ) 


( 195 ) 


( 196 ) 

( 197 ) 


By the same procedures as were used for deriving 
the velocity transfer, we find the components 


0 X = K* ((VK) 2 ) 

for the direct cascade and for the reverse cascade, re- 
spectively. 

The inertia turbulence for the direct cascade is 
governed by the spectral flow 


( 198 ) 

( 199 ) 
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I 


© x K' <^(l7.N , ) 2 / >= £ n , (200) 


where K' has been derived in (152) and (183) for the weak 
coupling and the strong coupling, respectively, and 
= ej (k =■ -) . 

We find the following results for the spectrum 
F Nv (k) of ^-fluctuations as follows: 


F N y (k) = 1.634 ( £ N y £ v 1/3 ) k~ 5/3 , zero coupling 


( 201 ) 


^ ^ I 

F N,v (k) = /=■(—) N ty / stron 9 coupling 


( 202 ) 


The parameter is 


N C,v- (£ N,AP 


1/2 


( 203 ) 


By an interchange of the subscripts v, u, we find 
the spectrum F N y (k) of N u - fluctuations as follows: 

F N , U 00 = 1*634 ( ^ )k” 5//3 , zero coupling * 204 ^ 


F N,u (k) 


= -J-(v/u f ) 2 N 2 k” 1 , strong coupling, 
J 3 * 1 ( ’ U 


(205) 


with 


N /,u = 


1/2 


( 206 ) 
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12. Summary and discussion 
12.1 Criteria of chaos 

A nonlinear dissipative system can be represented 
by the equation of evolution in the form (real or com- 
plex^ as follows: 

A C A 

O t + M) Y (t,x) = X(t,x) , (207) 


with 






(208) 


A 

X = 


X 


self 


S' 
+ X 


ext 


(209) 


where ? is an unknown random function (scalar or vector) , 

A 

is a non-random linear differential operator, and M njl 
is a differential operator that causes nonlinearity by 
involving ¥ or another random function, linear or non- 
linear. The driving force consists of the component X self 
that is reducible to 1 by a self-consistency condition, 
and the component X ext that is externally given and is 
random. The hydrodynamic equations of two phases, the non- 
linear propagation of light and of solitons, and the dy- 
namics of neurons fall into this category. 

The hydrodynamic equation is inhomogeneous and there 
is an advantage in transforming it into the homogeneous 
equation 
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(210) 


O fc + L)f (t,x,]ft = 0 


of higher dimensions, called the master equation. The 
differential operator is 

A A /*> \ 

L » M + X 0 , (211) 

with 3 = 3/3i|>. The random differential operator M in (207) 
is transformed into a non-random differential operator M 
in (210) . The field X that was a driving force in (207) 
becomes the only random convection in (210) . 

By writing 


f(t,x,t) = ^ (V - f(t,x)"] (212) 

and by taking the moment of (210) , we recover the hydro- 
dynamic equation (207) and show the self-consistency of 
the master equation. 

The hydrodynamic equation and the master equation 
can serve to describe the microdynamical state of tur- 
bulence. From the theory of universality of chaos, it 
will be interesting to know whether such a system can 
develop chaos and yield a broad spectrum. The theory of 
chaos is not yet able to predict the criteria. The hydro- 
dynamical closure by the degeneration of high orders into 
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lower orders can not find physical support. With the lack 
of a small parameter in turbulence, the kinetic method by 
Bogoliubov can not justify its closure. In our group- 
scaling method, we scale the master equation into a 
hierarchy of equations for the macro-group, micro-group, 
and sub-microgroup, representing the spectral evolution, 
the transport property, and the relaxation. The last 
process is represented by a network of random higher- 
order distributions. By the successive collisions, the 
distributions in the network lose their memory, and as 
an ensemble can be simulated by a porous medium which 
then closes the hierarchy by a self-consistent condition. 
The relaxation thus formulated makes the transport 
property approach its equilibrium. Hence, the kinetic 
equation is obtained. By the return to the continuum, 
we calculate the transport coefficients, and determine 
the spectral structure. 

We find a direct cascade and a reverse cascade. If 
one or both cascades can find their spectral balance 
from the equation of spectral flow, the chaos can de- 
velop into a broad spectrum. 

12.2 Dusty gas 

For a dilute suspension, where the particles are 
not sufficiently dense to produce a pressure, such that 
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«V u » E v 


(213) 


we simplify the system (19) -(21) for the v-phase into 


v A A A ^ A 

o + v V + U ) V = * u 

^ *** <*» V V ** 


(^> t + v-V)N v = 




(214) 


(215) 


with 


A 

< 


A 

= * N 
V V 


(216) 


while we keep the system for the u-phase in the original 
form (23) -(26). These two new systems consitute the 
'dusty gas' model by Saffman (1962). 

The spectral results for the dusty gas can be de- 
rived by a degeneration of our general formulation, under 
the condition (213). 

When we further neglect the coupling in (214) , we 
obtain the Burgers equation 


( 'b + V’Y - V V ) v = 


(217) 
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and the diffusion equation 


O t + v-V - = 0 . (218) 

The treatment of turbulence on the basis of the homogene- 
ous Burgers equation can be made by following the pro- 
cedures related to the master equation. 

12.3 Reduction of turbulence by the particles 
The power law 


K' = a” 1 I V u 1 k" 2 (219) 

u f I v - 

for the eddy viscosity has been suggested by Prandtl for 
turbulent boundary layers with a mean velocity gradient 
| Vu | . With this eddy viscosity and the equation of spec- 
tral balance in the form 


K' R 
u u 


- 


we find the spectrum 


F (k) = A u 2 k" 1 
u r r 


( 220 ) 


( 221 ) 


where u^, is 
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(£ u ''|?3I ) 


t 


( 222 ) 



1/2 


and is a numerical coefficient. The power law (209) 
confirms the earlier prediction by Tchen (1953, 1954). 
This law has been observed in the atmospheric boundary 
layers (Tchen et al., 1985) and in the tropopause (otten 
et al. , 1982) . 

It is interesting to note the analogy between (185) 

and (221) that can suggest that a shear flow with the 

combined effects of |Vu| and ct will have a reduced spec— 

— 1 

trum. This reduction can be seen from the lidar measure- 
ments in the atmosphere before and after the rain 
(Gurvich and Pokasov, 1972) . 


12.4 Equality of transport coefficients 

We choose a strong coupling where the governing 
parameters are a v , a u . By the use of Prandtl's formula, 
we write — 


K' 


v| v 



(223) 


and obtain 


V k) - * \ 2 *" 3 


( 224 ) 
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and similarly 


V k > - \ <** k * 3 • (225) 

The numerical coefficients are not determined in our di- 
mensional reasonings. These spectral distributions occur 
in a range where the transfer balances the production and 
where a v , a u are the only dominant parameters. 

We rewrite the formula (144) for the eddy viscosity 
in the form 


K' 


v| v 


and similarly 


K' 


u\v 



(k-)% 

F v (k") 


(k") 

9 * 


(k») 


/ 



(226) 


(227) 


Use of (166) is made for strong coupling. The subscript 
( ) u j y denotes that the transfer coefficient is endowed 
by the u' -fluctuations along the v-tra jectory , and ( ) v ^ v 
denotes that the transport coefficient is endowed by the 
v' -fluctuations along the v-tra jectory . Upon substituting 
for (224) and (225), we obtain 
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9 


(228) 


K ’v|v “ K ’ujv 


confirming the equality (8) from the Langevin method. 
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GROUP-KINETIC THEORY OF DISPERSIVE SOLITON TURBULENCE 
C . M . Tchen 

Department of Mechanical Engineering, City College, 

New York, N.Y. 10031 

As the fluid analogue, the nonlinear Schrodinger equation 
has a driving force in the form of emission of solitons by 
velocity fluctuations, and is used to describe the micro- 
dynamical state of turbulence. In order for the emission 

together with the modulation to participate in the transport 
processes (spectral evolution, eddy collision, and relaxa- 
tion), the non-homogeneous Schrodinger equation is trans- 
formed into a homogeneous master equation. By group- 
scaling, the master equation is decomposed into a system of 
transport equations, replacing the Bogoliubov system of 
equations of many-parti cl e distributions. It is in the 
relaxation that the memory is lost when the ensemble of 
higher-order distributions is simulated by an effective 
porous medium. The closure is thus found. The kinetic 

equation is derived and is transformed into the equation of 
spectral flow. The emission deposits the energy into the 

spectrum for accumulation toward small wave numbers by a 
reverse cascade with a negative eddy viscosity. The modula- 
tion develops a direct cascade and the modulation collision 
provides a fluctuation-dissipation for the drain. The spec- 
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tral laws for the soliton energy are k -1 , k"^ , and k" 6 in 

the inertia, coupling, and dispersion subranges, respective- 
ly. The spectral flow determines the statistical criteria 
for a nonlinear stochastic system to become chaotic. It is 
found that the homogeneous Schrodinger equation cannot de- 
velop strong turbulence. 

I. INTRODUCTION 

Many nonlinear phenomena in the continuous media exhibit 
fast density fluctuations superposed with slow ones, and the fast 
components have envelopes that fluctuate slowly. Such a modeling 
is called the soliton formalism and is usually represented by a 
homogeneous partial differential equation of the parabolic type, 
called the homogeneous Schrodinger equation. * This equation is 
considered more appealing than the Navi er-Stokes equation and has 
been widely used in plasmas, optics, and fluids, including atmos- 
pheric waves, gravity waves and sea surface waves. Unfortu- 
nately, the numerical computations of the homogeneous Schrodinger 
equation have not found chaos and deemed it not suitable for 
turbulence. From the hydrodynami cal standpoint, the Navier- 
Stokes equations for density and velocity can be transformed into 
a wave equation for the propagation of density with a coupling to 
velocity. Such a coupling is known as the production of noises 
by turbulence or the emission of solitons by velocity fluctua- 
tions, and must be represented in the Schrodinger equation to 
make it self-consistent with the Navi er-Stokes equations. The 
fluid analogue of solitons, i.e. the derivation of the 
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Schrbdinger equation from the Navi er-Stokes equations, has been 
treated by several authors. 2 ' 4 It is found that for the 
Schrbdinger equation to develop turbulence, the emission must be 
present. Hence, it is not justifiable to conclude that solitons 
cannot coalesce or break by collision on the basis of the 
Schrbdinger equation truncated to its homogeneous form. The 
dispersion in the differential form v 2 is important, without 
which the partial differential equation of Schrbdinger would be 
reduced into an ordinary differential equation in the class of 
Duffing oscillator. 

As the scope of treatment, we consider the non-homogeneous 
Schrbdinger equation to describe the microdynamic state of soli- 
ton turbulence, and transform it into a homogeneous master equa- 
tion (Section II). By a scaling into groups and by the loss of 
memory in the relaxation, we find the closure, and derive the 
kinetic equation of turbulence (Section 1 1 1 - V ) . By taking the 
moment we obtain the equation of spectral flow. By the emission, 
the energy is deposited into the spectrum, to be accumulated to- 
ward small wavenumbers by reverse cascade and transferred toward 
large wavenumbers by direct cascade. A transport theory deter- 
mines the eddy viscosity (Sections VI, VII). The equation of 
spectral evolution is solved (Sections VIII, IX). The group- 
kinetic consideration of spectral flow determines the criteria 
for a nonlinear dynamical system to develop chaos and turbulence 
(Section X). 
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II. KINETIC DESCRIPTION OF THE MICRODYNAMICAl STATE OF SOLITON 
TURBULENCE 

The nonlinear Schrbdinger equation 

(i + v> V 2 - i«AN)E = X (1) 

for the soliton field E(t,x) under the driving force X(t,x) has a 

nonlinear modulation by the density fluctuation 

N = - * A | E| 2 (2) 

2 

and a dispersion by V n V > where V> n » 

<0 and are real constants, having the dimensions of 

n 

viscosity, frequency and the inverse square of acceleration, res- 
pectively. The driving force is real, to be specified later. 

The random field 

E = E + E (3) 

v* 

can be decomposed into an ensemble average = <E> and a 

fluctuation E. The operators 

A = A + A (4) 

A 

can be used, where A = 1 is a unit operator, and 

A = 1 - A (5) 

is the deviation from the average. 

In the following, we assume homogeneous and isotropic turbu- 
lence in three dimensions, without mean field, i.e. 

1 = 0 . ( 6 ) 

By the distribution function 

f (t,x,E) = S [E - E ( t , x ) ] , (7) 

W, ^ 

the Schrbdinger equation in the physical space t,x is transformed 
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into the master equation 


O 

II 

<*♦- 

+ 

/O 

(8) 

in the phase space t,x,E. Here E is an 

'** ** 

independent variable, as 

A 

distinct from the function E(t,x) and 


A -V /v 

L = L v> + L N + L X 

(9a) 

L V S - 1 ^ n ** 

(9b) 

II 

(9c) 

L y « -1 X • * 

A w* w 

( 9 d ) 

are called differential operators, with 


\ = / > t , > * ^ A E . 

(10) 


It is not difficult to verify that the moment of (8) reproduces 


(1). The kinetic description (8) has the definite advantage over 
the fluid description (1) by being homogeneous. It is also the 
only logical way of developing a transport theory of spectrum- 
dependent eddy viscosity. 

The spectrum F^(k) gives the field intensity 

2^'Tk' F E (k') =<|£|'> . (11) 


It is defined by the Fourier transformation 

F E ( k ) [ d| e ■ i 3*?<|E‘(t,x)-E(x-i)|> 

L (27T) J - “ - - - (12) 

of the correlation function 

<Pt,x 1 ).£;(t,x 2 )|>= ]l d E j d E p E r E 2 <Rt l x 1 ,E 1 ) f(t,x 2 ,E )(> 

s JJ ^ E 2 d E 2 E ^-E 2 1 3 ) 

through the intermediary of the pair-distribution function 

^ 12 (t ,Xi , E i ; t,x 2 ,E 2 ) = <jf # (t,Xi,Ei) f(t,x 2t C 2 p. (14) 
Thus the determination of this spectrum becomes a problem of the 
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many-body statistical mechanics, involving a hierarchy of distri- 
bution functions. ^ 

Alternatively, by considering a finite upper limit of inte- 
gration k in (11), we get the spectral intensity 

d k * F E ( k * ) = < | E ° 1 2 > - (15) 

from which we find the spectrum F^(k) by a differentiation with 
respect to k. This method deals with the singlet -distribution 
function, or the moment 

E°(t,x) = f dE E f°(t,x E) , (16) 

or the energy 

<|E°| 2 > = f dE E -<|e*U,x) f° ( t , x , E )j> , (17) 

without involving the pai r- di s t r i b ut i on as was with the Bogoliu- 
bov method in (13) and (14). 5 

9 

The macro-group E is a component of the decomposition 

1 = E° + E ' . (18a) 


The micro-group 

E ' = E 1 + E " (18b) 

V* 

is re-decomposed into a first-order group £l and a subgroup 

E " = E (2) + E (3) + ... (18c) 

which comprises all the remaining higher-order groups. Likewise, 
we decompose the distribution into: 


f = f + f 
o 


f = f + f' 
f’= f 1 + f" 

f" * f (2) + / (3) + ... . 


(19a) 

(19b) 

(19c) 

(19d) 
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The operators 

A ° , A', A" (20) 

can be used for scaling into the macro-, micro-, and sub-groups, 
respecti vely . 

By the nonlinearity of the master equation, will be 

coupled to f 1 , and subsequently to f", in a sequence of equations 
of transport. Hence the present group-kinetic theory deals with 
the many groups 

f # , f 1 , f" (21) 

of the s i ng 1 e t - di s t r i b ut i on function, while the Bogoliubov theory 
deals with the distribution functions of many points: 3 
T*( t , x . , E ) (22a) 

^12 (t ’ x l’ E l‘ t » x 2‘ E 2 ) = <fU,x 1 ,E 1 )f(t,x 2 ,E 2 |)> (22b) 

• 

The three groups (21) represent the three transport proces- 
ses of spectral evolution, eddy viscosity, and relaxation. The 
closure is not based on the decreasing magnitude of the groups 
but on their decreasing coherence with correlation times 

r E > t E > t E» (23) 

such that the memory-loss in the relaxation makes the eddy coef- 
ficient (e.g. eddy diffusivity) approach its equilibrium and be- 
come adi abati c . This adiabaticity is not present in the high-order 
closure of distributions or of moments. 

III. DERIVATION OF THE KINETIC EQUATION OF SOLITION TURBULENCE 
By applying the scaling operators A, A # , A*, A^ 2 \ we 
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transform the master equations into the following equations of 
groups : 

( \+ A L)f = - A L* f* (24a ) 

O t + A L)f'= - l°f - A° L 1 f 1 (24b) 

C^ t + L)f* = - L*(f + f° ) - A^L^f^ • (24c) 

An operator (e.g. scaling operator, differential operator, col- 
lision operator, and evolution operator) applies to all functions 
which follow. 

In the right hand side, we find the collisions from the 
fluctuations of the higher-order nearest-neighbor group: 


- A L*f°= - A T f = S{7} (25a) 

_ AVf 1 = - A° L ' f ' = £'jf*} (25b) 

- A 1 [_ ( 2 ) f ( 2 ) . aV'T' = % f 1 , (25c) 

as related to the adiabatic ( non -f 1 uc t u a t i n g ) collision coeffi- 

cients 

e, ■&', e - (26) 

to be determined by a transport theory in Section VI. This 

transforms (24a) and (24b) into the kinetic equations 

(> t + A H ) f = - (27) 

U t + A # L ) f * = - L tf f + £'{ f°} (28) 

- 6 


for the probability function f and the macro-distribution f, res- 
pectively. In the following we shall deal with the kinetic equa- 
tion (28) only. The governing collision coefficient and its 
relaxation are determined by the system of equations 

(* t + A 1 7 ) f 1 = - L 1 *? + f') + t>"f 1 (29a) 

and 
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(29b) 


(* t + A " 1 ) f " = - L " (7 + f e + f 1 ). 

The function within the brackets is operated by the collision 

coefficient in the role of an operator for the propagation of 
memory . 

IV. EQUATION OF SOLITON TURBULENCE 


By taking the moment of the kinetic equation (28), we derive 
the following equation of soliton turbulence: 

( ® t - i 7 + ii^ n AN)E = -iX + 0, 

to 

with the stress 

J ’- J d tE B'K} • 

A subsequent multiplication by E and an average give the equa 
tion of spectral evolution 

< 1 E ° 1 2 > - T # 


in homogeneous and isotropic turbulence, where 


W 


of 


-i[i<X * E > -i<X • E >] 


is the coupling function, and 


o ow o 

T = - <| E • J > 


I 


dE E . V €/|<|E ( t , x ) f (t- T ) |>J. 

is the transfer function. The Lagrangian correlation 


<| E°*(t,x)f° (t-T ) | > 

V%> M 

contains a memory. Since & is derived from 


- S + L x 


wi th 


L » - ‘i ^ * I** . 


(30) 


(31) 


(32) 


(33) 


(34) 


(35) 


(36) 


(37) 


by (9), the collision coefficient can be written into the two 
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components 



= M) ♦ M } . 

(38) 

as due to 

/ / 

L and L , respectively. Similarly, we have 

N A „ 



t," = 'C N ' + tx 

(39) 

From (31) 

and (34), it entails 
c o 0 



t = Jn + } x 

(40) 

with 

T * - t n + T °x> 

(41) 


°N ■ /*{ { f ‘\ 

w 

(42a) 


£ ■ j *n Gi f, l • 

(42b) 


V. LOSS OF MEMORY IN RELAXATION FOR CLOSURE 

The equation of evolution (29a) for f^ is controlled by the 
collision 

S^f 1 } , (43) 

and is transformed by moment into an equation of evolution for E* 
with the controlling stress 

J 1 = f dE E £" ( t , x ) ^f 1 ^ (44) 

in analogy with (31). An effective collision would be 

5" (t,x ) f 1 , (45) 

e 1 1 — 

giving 

J 1 - f dE E f 1 

• S" (t.x) E 1 ( t , x ) . (46) 

Mr ~ * v ~ 

The collision coefficient that operates on f 1 in (44) keeps 
a memory by its dependence on E. If we wish to lose the memory 
in the relaxation, it suffices to suppress the E-dependence, cor- 
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responding to the fluid simulation 

£"(t>x,E)j } ~= ^*(t,x,E-0) 

l ^ N (t,x) , (47) 

•v 

since 

L " ( t , x , E = 0) = L jj,( t , x ) (48) 

by definitions (9a) - (9d). Hence we approximate (44) by 

J 1 = jl 
-W 

~^jj(t,x)E^(t,x). (49) 

Upon comparing (46) with (49), we identify 

£ e ff(t,x) = Efc(t,x). (50) 

A fluid simulation can be conceived if the diverse distribu- 
tions in the subgroup { 1 9 d ) undergo a large number of sub- 

collisions, so that the subgroup as an ensemble behaves like a 
porous medium without individuality from the separate E-depend- 

ence. The porous medium offers a Darcy coefficient *&". 

N 


With this approximation, we reduce (29a) and (29b) into 

o 


( * t + A 1 ^ ) fl * - L^f + f ° ) 


and 


(51) 

(52) 

(53) 


<\ + A " L N )f" = - L" (f + f\ f 1 ) . 

The differential operators 

- . ■ 

L N and L govern the evolution operators^, u,u, as follows: 

a , ^ 

( ^ + A 1 £ ) ”^(t,tj) = o (54a) 

( \ + A " L h> u N< t » t l ) a 0 (54b) 

(*^+AL)U(t,t 1 )=0. (54c) 

We conclude that the sequence of transport equations (28), 
(51) and (52) for the distribution functions 
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( 55 ) 


f . f 1 . f 

forms a closed system. 


VI. COLLISION COEFFICIENTS 

A. Darcy damping island collision coefficient 

The equation (52) for relaxation can be integrated into 

f" = -A "/ dt U N (t,t-X) L”( t -T) (f + f°+ f 1 ) , < 56 ) 

by means of the evolution operators (54b). Upon multiplying by 

L" and scaling by A 1 , we find the flux in the form 
N t 

- Akjf -A 1 / dx L^(t,x) A" U N (t,t-X )L“(t-T) (f + f + f 1 ). 


(57a) 


or in the adiabatic form 

. t-»Q> 


r v i 

-A 1 L "f " - «lx <L''( t,x) A" XJ (t,t-T ) L " { t -T ) > f 1 

N N — N U 


= " f 

S t 


(57b) 


Hence we derive the adiabatic coefficients 

= / dx <L^ (t,x) A" U N (t,t-T ) LjJ (t-T )> 


= - ( i O) ) 


2 0 


and 


. co 


)“ = J dx <N" ( t , x ) A" U N (t,t-X ) U N (t-x )>. 


(58) 

(59) 


B. Collision coefficient g* 

H 


The procedure of calculating from (51) by restricting to 


L^ is analogous to that for 16^. We obtain: 
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I 


f t-»» 

&N = { ) dx <L N {t 'V A ‘ ^ N ( t , t -"C ) L^(t-X )> 

= -(* tO,, ) 2 D' 

n N 

°N = f dx <N'(t,x) A 1 ^ ( t , t - T ) N'(t-T )>. 


C. Collision coefficient g' . 

A 


The integration of (51) by restricting to gives 

fl = at ^ (t,t-T ) (f + f ) 

Upon multiplying by and scaling, we get 

-A® l l f l = kj dt d(t) A ! U(t,t-T ) L l (t-t)Jf + fV 

*#%<*> 4 (t-x;>V 

■ G;(f‘<t-T)^ . 

The adiabatic coefficients are found as follows: 

r eo 

C’{} =J d-c <LjJ(t) A 1 ( t , t-T ) L^t-T)>{}- 

‘ 1 • - 1 ) 

= J dT <X 1 ( t , x ) A'li(t,t-T) X’(t-T)> 


D ' 

-X 

The coefficients 

r;(t,x,E) and D ' ( t , x , E ) 
are operators and generate memory. 


VII. TRANSPORT THEORY 

o 


A. Direct cascade T 


N 


From (34) and (42a), the direct cascade is written as 
T n = - / dE E i % <|E.*(t,x)f° (t--r)|> 

« - *jj <|E°| 2 >. 


(60) 

(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 


(67) 
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The collision coefficient 
< 0 

is negative by (60), and is not an operator since it is indepen- 
dent of E. The direct cascade performs a transfer toward smaller 

scales. 

o 

B. Reverse cascade T x 

From (34) and (42b), we calculate the reverse cascade in 


the form 


T x * -/«£ E ( )l>) . 


(69) 


It contains two Lagrangian correlations: one is the auto-correl a- 


t i on 


<X’(t,x)A'lfl(t,t-*r)X'(t-T)|> 


(70a) 


of the micro-fluctuations X' by definition (65), and the other is 


the cross-correlation 


^(T) = <| E**(t,x)A 4 U(t,t-T )f°(t-T )|> 


(70b 


«* t 

between E (t,x) and f (t-T ) 


By Fourier transformation, the correlation (70a) and the 
diffusivity (65) can be written in the form 

<X ' ( t , x ) A*'U(t,t-x)X'(t-T)> S' <X 1 (t ,x ) lf(t,t-T )X* (t-T )> 

^ ** — 4 vw 

= fdk" /X <X 1 (k" ) X ' ( - k " ) > h(x, k " , E ) (71) 

J W ^ ^ H _ A* Vi 

O' * f dk” *(X' (k")X* ( - k " ) > / dr h(T,k",E), (72) 

£ X V ** *+ ** — Jq ^ 

by using the orbit function 


h(r,k\E) = U(t,t-f) 

A 


(73) 


that is equivalent to the average evolution operator in the 



Fourier form. The Fourier transformation is made with respect to 
position in a length interval 2L in three dimensions within which 
the function X' is quasi -homogeneous , resulting in a factor of 
truncati on 

(tt/L) 3 . 

The transformation with respect to time is not necessary, because 
the temporal fluctuations are slow as compared with the spatial 
fluctuations in strong turbulence. The transformation with res- 
pect to both time and position would be necessary in weak turbu- 
lence to establish a resonance. 

C. Subdynamics 

Note that the Lagrangian function can be written in the 

form 

U(t,t--r )X 1 (t-T ) = X'Ct-T, x(t-x)], (74) 

where x(t-T ) satisfies the subdynamics of the detailed 
trajectory, as governed by 

x(t.) = E(t.) , E(t ) = - iX(t ) , (75) 

_ 1 - 1 - 1 «wl 


with the initial conditions at time t 

x ( t ) = x, x(t) = 0, E(t) = E, (76) 

and with t^ = t-~V. By analysing the trajectory, we find that 
the orbit function 

h(x,k'',E) = h '’r , k " ,E )h. (t, k " )hn(*c , k" )h (t, k " ) (77) 

" v — " wJ ~ 

has the following components 
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hfr(r,k M f E) = exp (-iik "• E T 2 ) 


( t, K" ) = exp ik" 2 ^T 

h (T, k" ) = exp fc" T 

" W /v/ 


h (T, k” ) 


z. <exp ik". ^ (t)> = exp^.T 


(78a) 
(78b) 
(78c) 

' N - . ( 7 8 d ) 

They refer to the streaming by E for the memory propagation, the 
dispersion by V , the Darcy damping , and the Landau damping by 
path perturbati ons X( T ). The average in (78d) can be calculated 
by means of a probability of transition ^(t,$) for the random 
pathi travelled during a time interval T. The probability func- 
tion satisfies an equation of transition with a damping rate 



By definition (65) or (72), the diffusivity is energized by 
< X * X 1 > for a duration of correlation that is determined by the 
evolution operator in the detailed fluctuating form or its aver- 
aged approximate form, through which the memory can propagate and 
be stopped for the closure and for the adiabaticity of the dif- 
fusivity. With this regard, we have represented the transmission 
of memory through E by the orbit component ^ . Other factors 
(the fluctuating t£ , the E-dependent probability of transition- ) 
may cause the memory to continue further, but are suppressed for 
cl osure. 

D. Reverse cascade by negative viscosity 

In analogy with the auto-correlation (79a), we write the 
cross-correlation (70b) in the form: 

< f 1 °(x)=f dk 1 *<|Ej (k')f*<-k',E)|> h(T,k',E). (80) 

The orbit function ^,(T,k ',E) has the same structure as (77). 
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We collect the results (72) and (80) and transform (69) into 
oo 

T °=/dr/d£ E i f d k " X (x •( k" ) X '( -k ••)> f dk • H° . (x , k " , k ' , E ) , (81) 

Xq J — I ' r ~ S ~ j “ TS1 — V- 

wi th 

H* (T , k 1 , k " , E ) £ } [h(T,k" ,E)) | E a . (k 1 ) f* (-k\E)|> h (t , k • , E ) ] . 

rsi w w, ^ - s 1 -w ^ 

(82) 

The two orbit functions 

h(T,k",E) and h(T,k\E) 

govern the small and large scales, respectively. They compete 

for the role of memory propagation by E. We attribute this role 

<*■* 

to the large-scale orbit function , by writing 

"a * h (x, k ' , E ) = - 1 k'k'x 4 h(x,k',E) 

•+- H* W, ***" k M H 

* - i k'k'x 4 (83a) 

h (t , k " , E ) *2 h (t , k " , E = 0). (83b) 

When the memory propagation is terminated, other roles (disper- 
sion and dampings) are taken over by h(x, k", E = 0) in view of 
the condition k'<<k". These considerations simplify (82) into 
H r S i(’ r »k' *k",E) =- h(x,k " , E = 0) ik ' k ' x 4 <|E*(k’) (-k • , E)| > , (84) 

and subsequently (81) into 

T x = ■ J ? <x r ' ( i" ) x; (-k")> 

dk ' k ' k' in 4 h (x,k" E=0) 

- r s y 0 - ' ~ 

dE E. ^<|E.° (k*) f®(-k' ,E)|>. (85) 

After integration with respect to E and rearrangement, we obtain 

Tx = -trace J dk" ;* <X ' (k")X • (-k")>[ T x (k") ] 5 R^ 

= - k*rJ . 

Here the vorticity function is 

= f dk • k • 2 X <| E # (k • ) | 2 > 


*/ 

*J 
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(87) 


= 2/‘ k dk , k' 2 F E (k' ), 

Jo 

the relaxation time is T (k"), such that 

r- CD 

C T ( k " ) ] 5 = real*/ dr r 4 h(r , k 11 , E =0 ) , 

X J o ~ 

and the eddy viscosity is 

K' = trace f dk"^< X'(k") X'(-k" ) >[T Y (k" ) ] 5 

J V— v* vv A 

dk" 


( 88 ) 


■ ! I 
J 0 


F x (k") [T- X (k")] 5 . 


(89) 

(90) 
(91) 


The spectral functions are F E (k) and F x (k), such that 

< I E^° I 2 > = 2 J k dk ’ F E ( k ' ) 

<X° 2 > = 2 A k d k " F (k") . 

- Jo * 

The transfer at a negative transport coefficient -K' indi- 
cates a reverse cascade. 

The relaxation time T(k"), as defined by (88), is calculat- 

Z 

ed by (78c) and (78d), as follows: 


C T x ( k " ) 3 ^ = y dr t 4 exp ( - - £„) 


- e.r 5 . 


N ” N ; 

Upon substituting (92) into (89), we obtain the eddy vis- 


cosity 


K * =4 


f 


dk" F y (k") (-£," - fc N ) 


-5 


(93) 


H 


E. Coupling function 

The coupling function, as defined by (33), can be inter- 

0 0 

preted as the flux of transport of E by X - fluctuation. Upon 

*** 

integrating (30), we write 


i/ dT A*U(t,t-T)X(t-T). 
<Jr\ - 


E~ = -i/ dT A"U(t,t-T) X (t-T). (94) 

/ 0 o 

Upon multiplying by iX and averaging, we determine the coupling 


f unct i on 
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« 


/"*> 


W = 3D^ , 

= f dk“ V <X*(k ' )• X* (-k ' )> / dr h (<*, k 1 , E=0) , (95) 

J -" ~ ~ _ ^>/ 0 — — 

with 

D° = trace D* > (96) 

by (72). The evolution operator U absorbs the damping effect of 

_ $ 0 

from J in ( 30 ) . 

o 

The orbit function that belongs to the trajectory of t is 

h (t.kj ,E=0) = exp (ik ,2 V n + 5^ + ^)t (97) 

by (77). A integration and a substitution into (95) yield 

W = ^ kdk ' V k,) ( "^N '^N } f w ( k'/k v ). (98) < 98 > 

The function 


^ = Real [1 - i ( k* /k^) 2 1 
= [1 + ( k '/ k ^) 4 ]" 1 

regulates the dispersion at the cutoff wavenumber 
k v * C( -«N -£*>/•>„]*. 


(99) 


( 100 ) 


VIII. MECHANISMS OF SPECTRAL EVOLUTION 

The evolution of the spectrum is governed by the transfer 

o O 

function Ty for reverse cascade, the coupling function W , and 
the transfer function T* for direct cascade. Two types of 

nonlinearity enter: the modulation by N and the emission by X. 
The spectral flow, as described by (32), involves cascade 
coupling, source and drain. 

A. Modulation and emission 

It is to be recalled that the soliton formalism was derived 
by a fluid analogue from the Navi er-Stokes equation for 
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compressible fluid, by considering a density fluctuation of two 

scales. The slow density wave is related to modulation in the 

detailed form (2), or in the form of intensities 

<N' 2 > - o< 2 [<| E' | 2 >] 2 . (101) 

This relation is called the equation of state. The fast density 

wave obeys an equation of wave propagation with the emission 

caused by velocity fluctuations. By an envelope transformation, 

the wave equation is transformed into the parabolic equation (1), 

called the Schrodinger equation for driven solitons, with the 

driving force X, such that 6,7 
«*• 

V . X ^ A 77 : u u (102) 

The propagation of finite sound wave and the excitation of 
acoustic turbulence are based upon the sawtooth shock wave 

formation having the time ^ as an invariant. This time is 

obtained by the effective slope in the propagation of a 

fluctuating velocity in the x- space. With this invariant as a 
parameter and by dimensional considerations, it is seen the 
velocity fluctuations follow the spectral law 

F u (lO ~ k -3 (103) 

It follows, from (102), the spectral intensity 

<X 1 2 > = T " 6 k" 2 (104) 

- s 

for the driving force. For the determination of x g , see Ref. 6. 
B. Instability as a source 

At the larger-scale end of the spectrum, the reverse cascade 
accumulates the energy from the micro-fluctuation into the 
macro-fluctuations, to be ultimately balanced by the instability 
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of soliton waves as a source in the form 


* E = i ^ <|E| 2 >, (105) 

so that the governing equation is 



o 

= K ' R £ i n inertial subrange, (106) 

from (86). 


At present, we shall not examine the spectral composition of 
the source; however, it would be required for analyzing the 
wave-turbulence interaction in the non-universal range of the 
spectrum. 

C. FI uc t uat i on - di s s i pat i on as a drain 

A nonlinear random system can be dissipative. The 
dissipation can be prescribed by an external parameter, such as 
the molecular viscosity. A "fl uctuation- dissipation" can be 
induced by density fluctuations and serves as a drain in the 
spectral flow. Our soliton formalism finds a drain from the 
fl uctuation-dissipation. 

0 

The transfer function T^ would play the role of a direct 
cascade in classical fluid turbulence, transferring energy from 

0 VO ' 

E into and finally down to the molecular dissipation. With 

the absence of molecular dissipation in soliton turbulence, the 
transfer function becomes a drain itself, by writing 

(107) 

I 

to balance the larger-scale coupling function W . The governing 
equation for this coupling subrange is 

W* - T* * o, (108) 
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from (32). The approximation of the scaled collision coefficient 
into an unsealed drain simplifies (67) into 

T* S - tf N <|E'| 2 > (109) 

and (98) into 

W fl = ( - 2 ^ ) " 1 2j k dk 1 F x (k‘) ^(k'/k^ ). (110) 

Here is defined by (99), with 

k v h [(-2Jr N )/y n 3*, dll) 

from ( 100 ) . 

It is to be remarked that the structure of the collision 
coefficient can be calculated from (60) and (61) for the 

determination of the probability function f(t,x,E) by (27). 

H. 

Since now ^ serves as a parameter, its explicit structure is 
not needed. 

D. Spectral flow 

By the use of (107), we reduce the equation of evolution 
(32) into 

*E * - \ * “ - T N 

* 2 ( -2C N ) 5 <X' 2 > 2/ k dk'k' 2 F E (k • > 

^ - /q L 

+ (2^)-! 2j k dk' F x (k«) 

+ E N <| E # | 2 >, with < 0. (112) 

This integral equation will be solved in Section IX. 

IX. SPECTRAL STRUCTURE 
A. Inertia subrange 

The inertia subrange is governed by the constant reverse 
cascade as described by the following integral equation* 
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£ E = 2C" 2 V 5 < l' Z> 2 / kdk ' k ' 2 F(r(k ' ) , (113) 

from (112). The solution for F £ ( k ) is found upon dividing by 
<X' 2 > and differentiating with respect to k, in the form: 

F E ( k ) = i(-2^ N r s ) 5 r s £ E k _1 . (114) 

Use of (104) has been made. It follows from (101) and (114) the 
densi ty spectrum: 

F N (k) = * 2 <l£'l 2 > (- 2 C, N T s ) 5x s r E k_1, (U5) 

We can introduce the notations 

= (-2^ n t. 3 ) - * (H6) 

k L = C( " 2g> N T s )?r s 5 (117) 

to scale the magnitude of the driving force and the source <5 r . In 
these notations, (114) and (115) become 

F E (k) = i« E 4 kj'V 1 (118) 

F„(k) =K 2 <|E'| 2 >^ e 4 k-V 1 . (119) 

An integration with respect to k yields the intensity of field 
f 1 uctuations 

< 1 E 1 l 2 > = i k £ ‘ 2 A k / k <- >. d20) 

and <N' 2 ) is given by (101). 

B. Coupling subrange (k/k v <<| ) 

The coupling subrange follows the inertia subrange at larger 
wavenumbers. From (112), we write the spectral flow in the form: 
(-2l? N )- 1 2/ k dk • F x ( k ') f w (k'/k y ) + £ n 2^ k dk* F E ( k ' ) = 0 , (121) 

The approximation 

^(k/k,) ~ I. f or k/K^« | (122) 

will be made. Upon differentiating with respect to k and by the 

use of (104), we solve for F E (k) in the form 


1 
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(123a) 


F £ (k) = 4^ k" 3 
or 

< | E ' | 2 > = 4 O 4 k’ 2 , (12 3b) 

where 4)£ is given by (116). 

From (101) and (123b), we derive the density spectrum 

<N • 2 ) = (4* 4> V k‘ 4 (124a) 

E 

and the spectral intensity 

F N (k) = <8*«feV k' 5 (124b) 

C. Dispersion subrange (k/k^ >> 1) 

The function for the dispersive cutoff is given by (99) in 


the approximate form: 

Z'vjU'/k^) = (k'/k p ) 4 . for k'/k p >> l . (125) 

Upon substituting for (125), we solve (121) to find the 
spectra 1 

F E (k) = 2 (k y ^) 4 k' 7 (126a) 

F N (k) = | [*- (k p <c/ E ) 4 3 2 k' 13 (126b) 

and the intensities 

\ (k p°E )4 k " 6 (127a) 

(w 2 ) = [ \ * (W* ] 2 k " 12 • (127b) 


All the spectral results indicate that the spectra are 
increased by the driving force through the parameter d E and by 
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the instability as a source through the length parameter 



X. DISCUSSION ON THE CRITERIA OF CHAOS 

A nonlinear dynamical stochastic system can be represented 
by the equation 

( ^t + L, + AL N )E(t,x) = Y ( t , x ) . (128) 

Here E(t,x) is the unknown random function without average, i.e. 

«*a 

(6). The advection consists 

of a non-random linear part and a random part L N that causes 
nonlinearity by being related to E(t,x) or being another random 
function. The driving force Y(t,x) may be of an external origin 
or be self-consistent with E(t,x). All the functions may be real 
or imaginary, scalar or vector, on one or more dimensions. Many 
stochastic problems in hydrodynami cs and optics can be described 
by this stochastic system. For solitons, the nonlinear 
Schrodinger equation (1) is identical to (128), by writing 

\ = -i X (129) 

and by using (9b) and (9c) as differential operators. 

The Duffing oscillator 

mx + (<r - Ax 2 )x = Y(t), (130) 

la* I 

with constant coefficients m, or , and ^3 , is reducible to the 

cubic equation 


. . -V * v /v 

(1 _y . \ bAN ) E = X 

A** 

by an envelope transformation, where N is given by (2), and 
are constant coefficients. The Korteweg-de Vries equation 

O t + u rr + A — — \ + v )u = o 
>x 2)x^ 


(131) 

y, b 

(132) 
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can be considered as an alternative to the Schrodinger equation 
for long waves. Here ^ and V are constants. 

A large amount of numerical effort has been devoted to 
calculating the soliton models, and found that the simple ones 
that are analytically integrable cannot develop chaos. This has 
lead to the belief that an integrable nonlinear equation cannot 
develop chaos, as exemplified by the Schrodinger equation, the 
Korteweg-de Vries equation, and the Burgers equation in their 
restricted form. This conclusion seems to be direct, and calls 
for the need of a causal relation between chaos and existence of 
solution, or for a statistical reason why certain dynamical 
properties are necessary to produce cascades and maintain a broad 
spectrum. 

The statistical theory as developed here may clarify to some 
extent this problem. As the starting point, we consider the 
dynamical equation (128) to describe the microdynamic state of 
turbulence. The equation has the difficulty of being non- 

homogeneous, with the consequence that the driving force will 
cause a hierarchy of high-order moments and is hindered to 
participate in the transport processes as does the nonlinear 
operator. To lift this hindrance, we use a kinetic method by 
raising (128) to a high dimension, and transform it into the 
master equation (8). The nonlinear operator then becomes 

t - + L n + Ly, (133) 

with 

L Y = Y • ^ . (134) 
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The group-kinetic method derives the kinetic equation of 

turbulence and establishes the spectral flow in the form (112). 

0 

The governing transport functions are the coupling functions W , 

0 

the transfer function T^ for the reverse cascade, the transfer 
0 

function T^ for the direct cascade, the source and the sink . 

The notation (129) for the solitons is followed. The solitons 

0 0 0 

gain energy from the pumping X into E by an amount W . The 
energy deposited is accumulated toward the small wavenumber end 

o 

of the spectrum by the reverse cascade T and is transferred 
toward high wavenumbers across the spectrum by the direct cascade 
At this end of the spectrum, the transfer function 

vHisV s -Sn<I!V> <135) 

serves as a drain in such a way that a collision coefficient is 
found and plays the role of a "dissipation" induced by 
turbulence, called "f 1 uc t uat i on- di ss i pat i on . " The energy balance 
between the reverse cascade and the instability on the one hand, 

and between the coupling and the drain on the other hand, even when 

the sink is absent, determines the two subranges and provides a 
broad spectrum of turbulence. 

If the soliton dynamics is described by the homogeneous 

A/ 

Schrodinger equation, i.e. X = 0, the spectral flow is governed 
by the direct cascade 

> t <|E°| 2 > ■ - Tj , with t * n >0 . (136) 

The cascade can not find a dissipation to complete the transfer, 
except in a decaying system, i.e. 
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t < | E*| 2 > > 0. 

The statistical considerations above help to determine the 
criteria of chaos . 
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